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Abstract 



We consider proton-nucleus collisions at high energy in the Color Glass 
Condensate framework, and extract from the gluon production cross- 
section the probabilities of having a definite number of multiple scatter- 
ings in the nucleus. Various properties of the distribution of the number 
JL . of multiple scatterings are studied, and we conclude that events in which 

^ the momentum of a hard jet is compensated by many much softer parti- 

(~| , cles on the opposite side are very unlikely except for extreme values of the 

saturation momentum. In the same framework, we also investigate the 
possibility to estimate the impact parameter of a proton-nucleus collision, 
, from the measure of the multiplicity of its final state 

Preprint CERN-PH-TH/2006-126, SPhT-T06/075 



1 Introduction 

Hadronic collisions at high energy involve the interaction of partons that carry a 
very small fraction x of the longitudinal momentum of the incoming projectile. 
Since the occupation number for such states in the nucleon wave function can 
become quite large, one expects that the physics of parton saturation |1I2I3| 
plays an important role in such studies. This saturation generally has the effect 
of reducing the number of produced particles compared to what one would 
have predicted on the basis of a pQCD calculation with parton densities that 
depend on x according to the linear BFKL (Balitsky-Fadin-Kuraev-Lipatov) 
| 4I5| evolution equation. 
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The counterpart of such a large occupation number is that one can treat the 
small- a; partons by classical color fields instead of particles. To that effect, the 
McLerran-Venugopalan model |6I7I8) is a hybrid description, in which the small- 
X partons are described by classical fields, and where the large- a; partons - fast 
and therefore frozen by time dilation - are described as static color sources at 
the origin of the classical fields, in agreement with the fact that small- a; partons 
are radiated by bremsstrahlung from the large- a; ones. Originally, the MV model 
dealt with large nuclei, with a large number of high- a; partons (the number of 
valence quarks is 3 A if A is the atomic number of the nucleus). In the MV model, 
the large- a: color sources are described by a statistical distribution, which they 
argued could be taken to be a Gaussian for a large nucleus at moderately small 
X (see also [5] for a recent discussion of this point). 

Since these early days, this model has become an effective theory, the so- 
called "Color Glass Condensate" (CGC) |l()lllll2j . Since the separation be- 
tween large x and small x is arbitrary, no physical quantity should depend 
on it. This arbitrariness leads to a renormalization group equation, the so- 
called JIMWLK (Jalilian-Marian-Iancu-McLerran-Weigert-Leonidov-Kovner) 
equation jl3ll4ll5ll6ll7ll8ll9ll0llltT^ . that describes how the statistical dis- 
tribution of color sources changes as one moves the separation between large 
and small x. This functional evolution equation can also be expressed as an 
infinite hierarchy of evolution equations for correlators [201) and has a quite 
useful - and much simpler - large Nc mean-field approximation |21j . known as 
the Balitsky-Kovchegov equation. 

In the collision of two nuclei at high energy, gluon production is dominated by 
the classical field approximation, and calculating it requires to solve the classical 
Yang-Mills equations for two color sources - one for each projectile - moving 
at the speed of light in opposite directions. This problem has been studied 
numerically in |22"23'24"25"26' for the boost-invariant case, with extensions to 
include the rapidity dependence |27I28| . But in fact, for collisions involving one 
small projectile - like proton- nucleus collisions ~ one can assume that the color 
sources that describe this small projectile are weak and compute the relevant 
amplitude only at lowest order in this source. When this is allowed, it is possible 
to obtain analytical expressions for amplitudes and cross-sections. This was 
done in a number of approaches for single quark or gluon production 29 3 0I31L 
l32l33l34l35l3HI37l38j . as well as for quark-antiquark production ,39,40 41,421^ 
I44l45l46l47l48j (see 0^1 for a review). In this paper, we are going to limit our 
discussion to the case of single gluon production. 

One of the main features of the gluon production cross-section in proton- 
nucleus obtained in the CGC framework is that it includes all the multiple 
scatterings on the sources contained in the nucleus. In this paper, we discuss 
the distribution in the number of these scatterings. In particular, we study how 
the momentum of a high-pj^ final gluon is balanced by the recoiling momenta 
of the struck nuclear color sources. This question has practical applications in 
discussing whether one could observe a loss of back-to-back correlations at high 
p± - for instance events with a single high-momentum jet in the final state - in 
collisions between a proton and a saturated nucleus. Another application of our 



2 



study occurs when one tries to relate the multiplicity and the impact parameter 
of the collision. And of course, one may also try to characterize the nuclear 
partonic content at low x from the distribution of "debris" that are produced 
in the collision with the proton. 

Note that the manifestations of the Color Glass Condensate on the back-to- 
back correlations have already been investigated in various approaches j47l5U| . 
by looking at the angular correlations between pairs of hard particles. This 
azimuthal correlation has been measured for deuteron-gold collisions by the 
STAR collaboration at RHIC [ST] , which observed that the pattern of azimuthal 
correlations in d-Au collisions is very similar to that found in pp collisions. In 
particular, it has a marked peak in the correlation function at 180 degrees - 
indicating that jets come in pairs. In the present paper, we address a different 
question, which involves most of the same physics: we do not consider the angle 
of emission of the particles, but instead we keep track of the number of recoils 
above a certain threshold /c™'" given a "trigger" particle with momentum k±, 
and we calculate the probability for having a given number of such recoils. 

Our paper is organized as follows. In section|21 we recall the CGC formula for 
gluon production in proton-nucleus collisions, and we also discuss the Glauber 
interpretation of this formula. In section|31 we show how to calculate the prob- 
ability P„ of having n scatterings in which the recoil momentum is larger than 
a certain threshold fc™™, when the produced gluon has acquired the momentum 
k± in the nucleus. This is done by constructing a generating function for these 
probabilities. In sectional we present numerical results for this distribution in 
the MV model, as well as simple analytical calculations that explain the most 
salient features; and in section we compare them with what happens in a 
model that has shadowing and geometrical scaling. Finally, in section El we 
discuss the possibility of estimating the impact parameter of the collision from 
the multiplicity in the final state. 



2 Gluon production in proton-nucleus collisions 

Here, we use without rederiving it the formula for the gluon yield obtained in 
PS) (Eq. (107)). According to this formula, the number of gluons produced per 
unit of transverse momentum and per unit of rapidity reads: 

dNg 1 f (fk±_ 2 



J^klC{k^)^p{q^-k^), (1) 



(Fqj^dy mn^ql J (27r) 

where ipp is the proton non-integrated gluon distribution, that we won't need 
to specify further in the following. The function C{k±), introduced in [22], is 
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the Fourier transform of a correlator of Wilson lines: 



C{k^) ^ I d^x^e^^---- 1:^1^' {U\0)U{x^)) . (2) 

U {x±) is a Wilson line in the adjoint representation of SU{Nc), evaluated in the 
color field produced by the sources that describe the nucleus, and the brackets 
• • • ) denote an averaging over these color sources. Equation ^ is accurate to 
the lowest order in the density of color sources contained in the proton, and to 
all orders in the color density of the nucleus. Thus, a way to picture its content 
is to say that a gluon of the wave-function of the proton travels through the 
color field of the nucleus before being produced. 

At first sight, it looks like the process taken into account by Eq. Q is a 
2-^1 process, in which one gluon from the proton (with transverse momentum 
Qj^ — k±) merges with a gluon from the nucleus (with transverse momentum 
k±) in order to produce the final gluon of transverse momentum Qj^. One might 
therefore be tempted to conclude that the Color Glass Condensate predicts the 
production of monojets in proton-nucleus collisions. However, this conclusion is 
too simplistic. The first reason is of course that transverse momentum is con- 
served in the CGC framework. This means that if the final gluon has acquired 
a large momentum k± while going through the nucleus, this momentum must 
come from the color sources present in the nucleus. In other words, if one sums 
the recoil transverse momenta of the sources struck by the propagating gluon, 
they must add up to —k±. 




Figure 1: Two possible scenarios for the recoiling scattering centers in the pro- 
duction of a high-p^ particle. Left: the large p± of the produced particle is 
compensated by many semi-hard recoils. Right: one recoil absorbs almost all 
the p±. 

Therefore, the real issue in order to conclude about the possible existence of 
monojets is whether the recoil momentum is shared among many sources (each of 
them acquiring only a small momentum) , or on the contrary absorbed mostly by 

^The function C{k±) is thus related to the Fourier transform of the cross-section between 
a color dipole and the nucleus. Hence, in the case of a quark-antiquark dipole, one can use 
this connection in order to relate proton-nucleus collisions and Deep Inelastic Scattering on 
nuclei |35|. 
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a single source (see Figure ^ for a cartoon illustrating the two situations). If the 
first scenario holds, then indeed one would have a high-g^ jet whose momentum 
is balanced by many soft recoiling particles - an event topology that would be 
close to one's idea of a "monojet" . In the second scenario, one would have a pair 
of high-gj_ particles, with almost opposite transverse momenta, in agreement to 
what perturbative QCD would predict. 

This interpretation in terms of multiple scatterings is particularly transpar- 
ent in the case where the distribution of color sources in the nucleus has only 
Gaussian correlations. This is the case of the McLerran-Venugopalan model (in 
which case the Gaussian distribution is local) , and also of the asymptotic regime 
believed to be reached after evolution to large rapidities with the JIMWLK 
evolution equation (in which case it is a non-local Gaussian distribution) |53| . 
Indeed, for a Gaussian distribution of nuclear color sources, it is possible to 
rewrite the function C{k±) in a form that has an obvious Glauber interpreta- 
tion. Let us reproduce here the main result of the appendix C of [SEl- Following 
Eqs. (C.5-6) of this reference, we can rewrite the function C{k±) as follows 

21 Z„-l 

X (27r)25(fci_L + --- + kn±- kMki±) ■ ■ ■ (J{kn±) . (3) 

In this formula, p is the number of scattering centers per unit of volume of the 
nucleus (assumed to be uniform), L is the longitudinal size of the nucleus, <7{k±) 
is the differential cross-section of a gluon with a scattering center of the nucleus, 
and (Ttot is the integral of the latter over kj_. Finally, /Xq = pL is the density of 
scattering centers per unit of transverse area. 



3 Distribution of struck scattering centers 
3.1 Definition 

In Eq. (PJ, the index n is the number of collisions of the gluon coming from 
the proton while it travels through the nucleus, and the exponential in the 
prefactor serves to unitarize the overall sum. Note that the integral over kj_ of 
the function C{k±) is equal to one, which means that this function should be 
interpreted as the probability for the gluon to acquire the momentum k± while 
going through the nucleus. The term of order n in this formula is therefore 
the probability that the gluon be deflected by a transverse momentum k± and 
undergo exactly n scatterings. By dividing this term by C{k±), we obtain the 
conditional probability that a gluon that comes out with a momentum k± has 
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scattered n times: 



2 L L L 



«i z„_i 

X (27r)2j(fei_L + • • • + fc„_L - ki_)a{kn_) ■ ■ ■ cr(fc„_L) ■ (4) 

So far, we have been a bit sloppy regarding the infrared behavior of the inte- 
grals over the transverse momenta that appear in Eqs. © and (0J. However, in 
the MV model for instance, cr(fe^) behaves as k'J^ at small fcj^ and it is necessary 
to introduce an infrared cutoff A in order for the integrals to be finite. It is well 
known that, although each integral behave as A"^, a partial cancellation occurs 
with the prefactor exp(— /Zgcrtot) so that C(A;_l) is only logarithmically sensitive 
to this cutoff^. Physically, this cutoff emerges from color neutralization that 
occurs on distance scales of the order of the nucleon size. Therefore, one should 
take A« Aqcd- 

This cutoff is of course also necessary in order to define the probabilities P^, 
so that they should in fact be interpreted as probabilities to have n scatterings 
with a momentum transfer larger than A. In the case of the Pn's, we can even 
push this logic further by defining the probabilities to have n scatterings with a 
momentum transfer larger than a certain fc™'" which is not necessarily related to 
A, and an arbitrary number of scatterings with a momentum transfer between 
A and fc™'". By doing so, we can explore how the distribution of the number 
of scatterings evolves with their "hardness". Let us denote Pn{ki_\k^^^) this 
probability. It is very easy to extract the relevant piece from Glauber formula, 
Eq. ©: 

p„(fcxifcr") = E j dz^j dz2--- j dz, 



p=0 



p+n 

Z\ 2p+Tl-l 



'^^ d'^ki± d'^kp± f d^fep+ix d'^kp+„± 



\ (2^)2 (27r)2 (27r)2 (2^) 

X (27r)^(5(fci^ H h kp+n±_ - k^)a{kij_) ■ ■ ■ cr(fep+„i) . (5) 

In this formula, n is the number of scatterings with momentum transfer larger 
than fc™'" and p the number of scatterings with momentum transfer between A 
and fc^''". 

■^In the individual probabilities P„ however, this cancellation does not occur and one has 
a quadratic sensitivity to the cutoff A. 
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3.2 Generating function 

Although a direct numerical evaluation of Eq. lO is in principle feasible, it turns 
out to be easier to compute the following generating function instead: 

+ OC 

F{z, k^ikrn = Pnik±\kTn 2" ■ (6) 

From this function, it is straightforward to go back to the probabilities P„ by 
the following formula:"^ 

Pn{k±Wr) - / ^ e-"'' F{e'', k^\ktn ■ (7) 



Therefore, it will be sufficient to calculate the generating function for complex 
z's on the unit circle. In practice, one should evaluate the generating function 
for a finite number (usually a power of two) of values z = e** , with the angles 
6 equally spaced on the circle, and then evaluate the Fourier sum by the fast 
Fourier-transform algorithm. 

It is easy to replace P„ by its expression in Eq. ||HJl, and to perform the sum 
explicitly. In order to disentangle the various variables fc^^, one must replace 
the delta function by its Fourier representation. This leads to: 



F{z,k^\k'r) ^ J <fx^e-^^---- 



X exp < ^0 



A fcfin 



(8) 



Note that, for z = 1, the numerator of this formula is identical to C(fcj_). This 
was of course expected, since fc^jfc""") = 1 (because this is the sum of all 
the probabilities P„). As one can see, the only difference between the calculation 
of C{k±_) and of the numerator in Eq. © is that the exponential exp(iZ_L • x±_) 
is weighted by a factor z for the values of l±_ above fc™'". Therefore, calculating 
the generating function can be done via a fairly minor modification^ of the 
numerical methods used in order to calculate C(fcj^). 

^Another approach to obtain the probabihties from the generating function is to compute 
the successive derivatives of the generating function at 2: = 0. However, this would require to 
evaluate derivatives of high order, which is very difficult to do numerically. 

■'This observation also indicates how to construct the generating function for probabilities 
that are more general than the ones considered here: in order to compute the probabilities 
Pn{n) to produce n particles in some region fJ of the single particle phase-space, one must 
weight the exponential exp(iZx -asx) by a factor z when Zj^ 6 f2. This approach could be used 
in order to study the recoils in a specific angular sector for instance. 
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In fact, as one can readily see, in order to calculate the argument of the 
exponential in Eq. (jSj), it is sufficient to compute the following two integrals, 

A 

B{x^\kT'')^ I ^^''"■'"^(^^) ' (9) 

as a function of x± and fc™™. 

In actual numerical calculations, the lower limits, at = A in and at 
l± = fc™™ in B, are implemented by multiplying the integrand respectively by 
'd{l±/A) and i?(/j_/fc™™). The function d{x) interpolates between at small x 
and 1 at large x, the transition between the two regimes being located around 
X = 1. One could in principle take for -dlx) the ordinary step function 9{x), 
which corresponds to sharp lower limits - as written in eqs. (jSjl - but such a 
choice generally leads to an oscillatory behavior of the functions A{x±) and 
i?(a;_L|fc™™) as a function of xj_. Choosing a function ^{x) that has a smooth 
transition between and 1 is helpful in order to tame these oscillations. 

Once the integrals A and B have been calculated, one can write: 

C(fci) = J (fixx_ e-'''^-''^ et^lM^^) ^ (iq) 

and then 



3.3 Models for a{k^) 

In the rest of this paper, we consider two different models for the differential 
cross-section a{l±). 

The first of these two models is the McLerran-Venugopalan model |6I7I8| . 
which assumes a local Gaussian distribution of color charges in the transverse 
plane for the nucleus. It is well known that this leads to^ 



o{l,_) = . (12) 



In the MV model, one can have important rescattering effects (tuned via the 
density parameter /ig), but there is no leading- twist shadowing. Note that the 
saturation momentum Qa is given by: 

\ QCD I 



^Here, the formula has been written in the adjoint representation, since it is a gluon that 
propagates through the nucleus. For a quark, one would simply have to replace the color 
factor by Cf = (A^^ - l)/2Wc. 
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Here, we have written the saturation momentum in the fundamental representa- 
tion, in order to facihtate the comparison with the values of Qa extracted from 
Deep Inelastic Scattering at HERA. 

The second model we will consider is based on a Gaussian effective theory 
that describes the gluonic content of a nucleus evolved to very small values of 
X, discussed in It corresponds to the choice 

In this model, hereafter referred to as the "asymptotic model", c « 4.84 and 7 
is an anomalous dimension whose value is 7 « 0.64. One of the peculiarities of 
this model is that it has the property of "geometrical scaling" , since it depends 
on the momentum and on Qs only via the ratio l±/Qs- Contrary to the MV 
model, this non-local Gaussian model has significant leading-twist shadowing, 
whose strength is controlled by the anomalous dimension 7 (more precisely by 
the departure of 7 from 1). 



4 Results in the MV model 
4.1 Multiplicity distribution 

Let us first start by displaying some results in the MV model. In Figure |21 
we first show the distribution of the probabilities P„ as a function of n, for 
Qa=2 GeV^ and various values of the threshold momentum fc™'". One can 
see that the width of the multiplicity distribution decreases with an increasing 
/c™™. This is of course quite natural, since by increasing A:""" it becomes less 
and less likely to have events in which there are a large number of recoils. Note 
also that for all fc™'" such that Qs <C fc™'" < k±, the most likely number of 
recoils is n = 1, while for k± < fc""" the most likely situation is n = 0. 

We can in fact understand analytically this distribution P„ in the situation 
where the momentum exchange kj_ between the incoming gluon and the nucleus 
is much larger than the other scales, Qs, fc™™ <C k±. This means that we need 
only to estimate the functions A and B defined in Eq. ^ for values of x± that 
are much smaller than the inverse saturation momentum, x± <C Qj^, and much 
smaller than (A:™'")"-'^ . This allows us to expand the exponential exp{il± ■ x±) in 
order to evaluate the integral over l±, leading to the following approximations: 

A{x±) ~ —- — - x\ In 



167r V^2;_L 



8tt 



1 a;^ . / 1 



- ^ In 



(15) 



2 vfcT'"2;_L, 

Then, in order to evaluate the generating function via Eqs. (|lUf) and Hll|l . one 
can use the following result, valid at large kj_, 

fx^ e-*'=-=^- e-^-- in(.V-^) ^ ^ . (16) 
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Figure 2: Distribution of the probabilities P„ in the MV model. The produced 
particle has acquired a transverse momentum kj_ = 10 GeV and the saturation 
momentum is set to the value Ql = 2 GeV^. The threshold fc™'" for counting 
the recoiling particles takes values between 0.4 and 16 GeV. The black solid line 
represents the result of an approximate calculation (see text). 



(The value of the constant has no influence on this result in the limit of large 
A:_L.) Thanks to this formula, we obtain immediately 

(z-1) 



i^(z,fcj^|fc7'^) « z e**"'*^™)"^ (17) 

One can see that in this limit, the generating function is universal in the sense 
that it does not depend on the momentum k± acquired by the incoming gluon. 
Moreover, the probability of having zero scatterings with a recoil above fc""", 
Pq = F{z = 0, A:_L|fc™'"), is zero. In other words, when /c""" ^ k±, there must 
be at least one scattering above fc™™ in order to give such a large kj_ to the 
incoming gluon. 

We can go a bit further, since it is easy to recognize that the generating 
function obtained in Eq. lfT7|l corresponds to the following distribution of prob- 
abilities: 

Po - 0, 

P„ ^ ^!4tt e-" with n = ^y^^l, . (18) 
(n-1)! 87r(fc™")2 ^ ^ 

In other words, the distribution of multiplicities is a Poisson distribution shifted 
by one unit. The physical meaning of this shift will become transparent later 
in the discussion. In Figure |3 we have compared for fc™'" = 0.4 GeV the 
numerically evaluated P„'s with such a shifted Poisson distribution, and as one 
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can see the two agree extremely well (except for Po , which is very small but not 
exactly zero). 

Note however that the value of n we had to use in this fit differs by about 
25% from the predicted value given in Eq. H18(l . This kind of deviation is 
expected, because this formula for n is only valid for fc™'" » A, a condition 
which is at best marginally satisfied for A;™'" = 0.4 GeV (we have taken the 
infrared cutoff to be A = 0.2 GeV). Moreover, from the approximations that 
have been used in order to obtain Eq. H15|l . a generating function of the form 
F{z) = z exp(n(z — 1)) ~ that leads to a shifted Poisson distribution - is 
obtained as long as k± ^ Qs,fc™™. It is only the accurate prediction of the 
value of n that requires in addition fc""" A. This explains why, despite the 
fact that n was not very accurately predicted at A:™'" = 0.4 GeV, the obtained 
distribution was nevertheless of the form given in Eq. H18() with a very good 
accuracy, because kj_ — 10 GeV. 

4.2 Number of recoils 

MV model - Qs" = 2 GeV^ 

10 r- 

1 r 

-S 0.1 r 

o 

OJ 

° 0.01 r 

B 

Z 0.001 r 

le-04 r 
le-05 - 

0.1 1 10 100 

K,"'" (GeV) 

Figure 3: Number of recoiling scattering centers in the MV model, as a function 
of the threshold fc™™. The saturation momentum is set to the value = 
2 GeV^, and the momentum k± of the produced particle is varied between 1 
and 50 GeV. 

Next, we display in FigureOthe average number of recoils above the thresh- 
old fc™™, defined as 

oo 

iV(fc^|fc';^'") = ^nF„(A;^|A;f") , (19) 

n=l 

as a function of fc""", for various momenta k± and a fixed = 2 GeV^. We see 
that the number of recoils grows significantly at small A""", and tends in this 




11 



region to become universal and independent of k± . Moreover, a striking feature 
of this number of recoils is that it is very close to unity for any value of fc™™ 
such that Qs <C fc™'" < k±. This means that when the gluon acquires a large 
momentum k± from the nucleus, there is always one hard recoil - and only one 
- that provides most of this large momentum. 

Again, it is possible to have an analytic understanding of these properties 
of the average number of recoils from Eq. H17(l . Indeed, the average number of 
recoils is given by the derivative of the generating function at z — 1, and we 
obtain 



Nik^lk'l"') « 1 + 



(20) 



87r(fc™'^)2 ■ 

This analytic expression is also displayed in Figure 13 and it reproduces well 
the numerical calculation for A;™'" ^ k±. It deviates from it at very small fc™™ 
due to a non-trivial interplay between fc™™ and the infrared cutoff A, which is 
not correctly captured by our simple analytic calculation. And of course this 
analytical result does not work for fc™'" > k± because this is outside the range 
of validity of our approximations. 



4.3 Momentum distribution of the recoils 



MY model - = 2 GeV" 



.a 

E 




0.01 



0.001 : 



le-04 



K,""" (GeV) 



Figure 4: Distribution of the number of recoiling scattering centers in the MV 
model, as a function of the threshold fc™™. The saturation momentum is set 
to the value = 2 GeV^, and the momentum k± of the produced particle is 
varied between 1 and 50 GeV. 

In Figure ^ we have taken the derivative of the average number of recoils 
with respect to /c™'", in order to obtain the momentum distribution of these 
recoils. At large k±, one can clearly see that this distribution consists of two 
components: a universal (almost independent of k±) semi-hard component made 
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of recoils with momenta of the order of Qs or smaller, and a component peaked 
around fc™'" — k± . The latter peak has an area unity, and it is simply translated 
when k± is changed. By taking the derivative of Eq. (|20|l . one can readily obtain 
a contribution that reproduces well the numerical result in the semi-hard region 

dN _ g^N^f^l ^2^^ 



dln(fcf'") 47r(fc™")2 ' 

In fact, it turns out that it is also possible to estimate this derivative in the 
region where /c""" is comparable to or larger than k± (both of them being very 
large compared to Qs). When both fc""" and k± are large compared to Qs (i.e. 
to /ip), it is enough to expand the exponentials of A and B in Eq. to first 
order, and write^ 

/ cPxi_ e-^'^^-'^-L \a{xi_) + {z- l)B(x_L|fc2^''') 
Fiz, fcilfcf ") « ^— ^ ^ ^ . (22) 

The multiplicity being the derivative of _F at z = 1, we have 

Going back to the form of A and B, we see that the integration over x± 
simply produces a S{k± — l±), making the integral over l± trivial as well. In 
this kinematical region, we obtain an extremely simple result: 

Nik^lkfn « = ^(fc^/fcf ") . (24) 

We see that this component of the multiplicity is nothing but the cutoff function 
that we are using in order to separate the momenta that are below fc™'" from 
those that are above. Therefore, the precise shape of the average number of 
scatterings for fc™™ above k± is not a property of QCD, but merely reflects 
the fact that we have an extended rather than a sharp cutoff. Nevertheless, the 
interpretation of this contribution is quite straightforward: when fc™'" is smaller 
than k± there is one recoil (that absorbs most of the momentum kj_), but it 
is unlikely that there is a recoil with a momentum bigger than the momentum 
k±. Taking a derivative with respect to fc™'", we obtain the corresponding 
contribution to the momentum distribution of the recoils: 

« = ^'(fcx/fcr) . (25) 



<iln(fc'["") k 

In Figure^ we have represented for k± — 50 GeV the sum of the contributions 
given in Eqs. H21|l and (|25|l (taking for the latter the same "step function" i9(a;) 

^Note that the 1 in the Taylor expansion of the exponential does not contribute at large 
k± since it only gives a term proportional to S{k±). 
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as the one used in the numerical calculation of the integral B). The sum of 
these two components reproduces with a fairly good accuracy the numerical 
results for all fc™'" down to fc™'" ~ 500 MeV. The small discrepancy between 
our analytical estimate of the peaked contribution and its numerical value is 
due to rescattering corrections - indeed, our derivation of Eq. retains only 
the leading-twist contribution. As one can see, the numerically obtained peak 
is slightly shifted to the left of the analytical result. This is easy to understand: 
since there are a few semi-hard scatterings in addition to the hard one, the hard 
scattering needs to provide a little less than the momentum k± acquired by the 
gluon. This shift is a form of coUisional energy loss (for a cold nuclear medium). 

Note that, when the "step function" ■d{x) becomes a real step function 9{x), 
Eq. (^5)1 would imply a peak proportional to S{k™™ — k±). However, we expect 
that higher- twist corrections to Eq. (|22|l would be important in this limit, and 
they are likely to smear out slightly the delta peak. 

Before considering the Qs dependence, let us come back to the Poisson dis- 
tribution shifted by one unit found in Eq. (|18f) . when fc™'" ^ k±. The shift by 
one unit is due to the fact that, when the threshold fc™'" is so low compared 
to k±, there is always at least one scattering (moreover, we know now that this 
scattering has a recoil momentum which is close to A;_l). The meaning of Eq. p8|l 
is therefore that the remaining n — 1 - semi-hard - scatterings that come along 
with this hard scattering have a Poissonian distribution, which merely reflects 
the fact that they are independent from one another. 

4.4 Dependence on Qs 

Finally, let us have a look at the dependence on the saturation momentum. For 
this, we set the momentum k± acquired by the gluon to 10 GeV, and we study 
the momentum distribution of the recoils for various values of Q^. The results 
of this analysis are displayed in Figure El As long as the saturation scale Qs 
remains small compared to k±, only the semi-hard part of the distribution is 
affected by changes of Qs , while the peak around fc""" = k± remains unchanged. 
The latter result is due to the fact that, since this peak is well approximated by a 
leading-twist calculation, it must be independent of saturation physics with the 
same accuracy. It is only when Qs becomes very large that one cannot neglect 
higher- twist corrections at large fc^'", and that the peak at fc""" = k± eventually 
disappears. Since the distribution of semi-hard recoils is quite sensitive to the 
value of Qs (it is proportional to g^^Q, which is proportional to Q^ up to a 
logarithm), it could perhaps be used as a way to estimate Qs- 

The disappearance of the peak also provides a qualitative answer to our ini- 
tial question regarding the possible existence of monojets: any parton produced 
with a k± which is much larger than the saturation momentum in the nucleus 
must have its momentum balanced by another parton on the opposite side (the 
latter comes from the scattering center that has undergone the hard collision). 
But all the partons with a transverse momentum comparable to or smaller than 
Qs need not have their momentum balanced by a leading parton on the oppo- 
site side, since it can be balanced by several softer particles (coming from the 
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Figure 5: Qs dependence of the distribution of the number of recoihng scattering 
centers in the MV model. The saturation momentum is varied between = 
1 GeV^ and = 20 GcV^, and the momentum k±^ of the produced particle is 
set to 10 GeV. 

semi-hard component of the distribution of recoils) . As long as the saturation 
momentum remains relatively small, say Qs ~ 1 — 3 GeV, this conclusion is not 
going to alter one's common expectations regarding jets: all hard jets with a 
momentum larger than say 10 GeV must come in pairs. It is only for a very large 
Qs that one would start seeing non-conventional event topologies where a hard 
jet would have its momentum balanced by a large number of softer particles. 




5 Effect of leading-twist shadowing 

Let us now briefly compare the results previously obtained using the MV model, 
with those one obtains by using the model defined by Eq. (|14|) . Basically, the 
two models - at an identical Qs - differ by the nature of the correlations among 
the color charges in the nucleus. In particular, the MV model does not have any 
leading-twist shadowing, while the second model has an anomalous dimension 
7 different from unity and thus provides some shadowing. It is believed that 
the latter model is a better description of a nucleus at very small momentum 
fractions x. 

In Figure El we first compare the average number of recoils for the two 
models (thin lines: MV model - dots: asymptotic model). The value of the 
"trigger momentum" k±_ is held fixed at a value of 10 GeV, and the saturation 
momentum squared is varied in the range 1 — 10 GeV^. One sees that the 
number of semi-hard and soft recoils is quite smaller in the asymptotic model 
than in the MV model. At the largest of the considered Qs, the number of soft 
recoils is ten times smaller in the asymptotic model than in the MV model. We 
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Figure 6: Qs dependence of the number of recoiling scattering centers in the 
asymptotic model, compared to the MV model. The saturation momentum is 
varied between = 1 GeV^ and — 10 GeV^, and the momentum k± of the 
produced particle is set to 10 GeV. Lines: MV model. Dots: asymptotic model. 



interpret this as an effect of shadowing, which "hides" the scattering centers 
from the passing giuon. A similar observation was made in . where it was 
seen that the multiple scatterings that lead to the Cronin effect are almost 
inexistent in this asymptotic model. Also, an effect of shadowing is that the 
dependence on Qs is much weaker in the asymptotic model: piling up more 
and more color charges in the nucleus does not lead to many more scatterings 
if the gluon cannot see them because of shadowing. This weaker dependence 
on Qs is also seen at large fc""", where one can hardly see any change even at 
g2 = 10 GeV^. 

Another feature of the asymptotic model is that the "plateau" at = 
1 for Qs <C fc™'" < k± is no longer really flat. Instead of a wide plateau 
between Qs and k± , one has instead a slow but steady rise of the multiplicity as 
^mm (decreases. For this reason, we expect the two-component structure of the 
momentum distribution of the recoils to be less pronounced in the asymptotic 
model than in the MV model. This is what we check by taking a derivative with 
respect to A:™'", as illustrated in Figure In these plots, one can see that the 
dip between the low-momentum component and the peak around A:™'" = k± 
is not as deep as in the MV model. This means that one should expect the 
distribution of momenta in the "away-side jet" to be more extended towards 
softer momenta, as one probes the nucleus at smaller and smaller values of x. 
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Figure 7: Qs dependence of the distribution of the number of recoiling scattering 
centers in the asymptotic model, compared to the MV model. The saturation 
momentum is varied between = 1 GeV^ and — 10 GeV^, and the mo- 
mentum k± of the produced particle is set to 10 GeV. Lines: MV model. Dots: 
asymptotic model. 
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6 Measuring the impact parameter 
from the multipUcity? 

Based on the above study, one can address a related question:^ is there a cor- 
relation between the measured multiplicity in a pA collision (event by event) 
and the impact parameter of the collision? and with what accuracy coiild one 
determine the impact parameter based on this correlation? 

In this theoretical study, the question one can answer is the following: if the 
measured multiplicity in an event is n (in addition to the hard jet of momentum 
k±), what is the probability distribution of the various impact parameters? In 
order to answer this question, we will make three assumptions: 

(i) When two bunches of nuclei and protons collide in an accelerator, all the 

impact parameters b are equally probable. 

(ii) The only recorded events are those where |6| < R, where R is the radius 
of the nucleus (we neglect the radius of the proton). Assuming here for 
simplicity that the trigger efficiency is the same for all 6's, the probability 
of a given impact parameter (in the absence of any other information 
about the collision) is a priori equal to P{b) = I/ttR? . 

(iii) The density parameter at a given impact parameter b is proportional to 
the thickness of the nucleus at this impact parameter, i.e. to V-R^ — • 

Let us introduce the probability 'P{n, b) of having simultaneously the impact 

parameter b and the multiplicity n (it is implicit in all this section that we 
mean the multiplicity above a certain threshold fc™'" when the passing gluon 
has acquired the momentum k±_ - these variables will not be written anymore 
in order to avoid encumbering the notations). P{n,b) must be normalized so 
that one has 

oo ^ 

J2 d^bV{n,b) = l. (26) 

n=0'' 

The probabilities P„ defined earlier in this paper can be obtained from this more 
general object by 

^ nn,b) 

" E„nn,6)- ^''^ 

The denominator is necessary so that the P„'s add up to unity. Obviously, this 
denominator is a function that depends only on 6, whose integral over b is unity. 
It is nothing but the probability of having a collision with impact parameter b, 

^This question is reminiscent of tiic attempts to measure the impact parameter in colli- 
sions on nuclei by counting the so-called "gray tracks". Usually, in the relatively low-energy 
collisions where this has been used, the picture is that the passing projectile would kick nu- 
cleons out of the nucleus and that by counting these nucleons one could estimate the impact 
parameter. The general idea of our study is the same, except that the action takes place at 
the partonic level. 
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when the incoming gluon has been scattered off the nucleus with a momentum 
fc . It is easy to convince oneself that this probabihty is given by: 

where the b dependence of C{k±) comes impUcitly via the parameter If 
there were no trigger bias, this quantity would simply be uniform and equal to 
{ttR"^)^^. However, because it is slightly more likely to have a large fc_L in central 
collisions than in peripheral ones, the mere fact of selecting a specific k± in the 
final state introduces a certain bias in the distribution of impact parameters®. 
Therefore, one has 

and we see that no new calculation is necessary. It will be sufficient to calculate 
Pn at fixed n as a function of b (the b dependence comes in via fj,Q ^ V — b"^) ■ 
From this object 'P{n, b), it is easy to obtain the normalized distribution of 
impact parameters conditional to having an event with the multiplicity n, which 
is the solution to the question we asked, 

We have evaluated this quantity numerically in the MV model. The only extra 
parameters that need to be set are the coefficient of proportionality between 
the density Hq and the size \/R^ — b"^ - we set it so that the saturation scale 
at the center of the nucleus (6 = 0) is 2 GeV^ - and the nuclear radius, taken 
to be i? = 6 fm. The results are displayed in Figure |H1 for events where the 
gluon acquires the momentum fcj^ = 10 GeV and with a threshold momentum 
of fc™'" = 0.4 GeV for counting the number of recoils. 

The results are fairly intuitive: events with a low multiplicity are dominated 
by large impact parameters, and events with a high multiplicity are much more 
central. But we also see that selecting a given final multiplicity only gives a 
fairly vague idea of the impact parameter, since the distributions of probability 
for 6 at a fixed n arc quite wide, with important overlaps between the curves 
for different final multiplicities. And to make things even more difficult, the 
two extreme values of the final multiplicity [n — 1 and n = 16 in our example), 
which have the least overlap in &, correspond to very rare events as one can judge 
from the figure |21 Therefore, it seems realistic to make two centrality classes, 
reasonably well separated in impact parameter, based on the observed number 
of recoils. Changing the value of the threshold A:™'" may help this separation, 
but we have not investigated that approach systematically here. 

^One can check numerically that this bias is significant only for very peripheral collisions. 



19 



MV model - K, = 10 GeV , K,"'" = 0.4 GeV 
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Figure 8: Probability distribution of the impact parameter for various final 
multiplicities, in the MV model. The trigger momentum is set to k± = 10 GeV, 
the threshold momentum to fc""" = 0.4 GeV. 

7 Conclusions 

In this paper, we have calculated the distribution of the number of scatterings 
in proton-nucleus collisions, in the Color Glass Condensate framework. This 
has been done by calculating the generating function for the probabilities of 
having a definite number of scatterings. We observe that, when the produced 
gluon has a transverse momentum which is large compared to the saturation 
scale, then this momentum is provided mostly by a single scattering center in 
the nucleus, leading therefore to the familiar di-jet configuration. This hard 
scattering is accompanied by a larger number of semi-hard scatterings, with 
transferred momenta of the order of the saturation momentum or smaller. By 
comparing the McLerran-Venugopalan model with a model that describes the 
regime of very small x, we also see that the shadowing present in the latter tends 
to suppress these semi-hard scatterings, and to blur the separation between 
the hard and semi-hard scatterings. Finally, we have discussed the correlation 
between the final multiplicity and the impact parameter, and shown that it is 
not a very strong correlation, that can at best be used to make a gross separation 
in at most 2-3 centrality bins. 

As a final note, let us mention that the results discussed in this paper are a 
particular case of some general results on random walks (in two dimensions in 
our case) where at each step one may have a random step size (both in magni- 
tude and direction), according to a certain probability law. If this probability 
distribution for the step sizes is falling very quickly, then the only way that the 
random walk may end far away from the origin is to add up a very large number 
of small steps. On the contrary, if this probability distribution has an extended 
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tail at large step sizes, such that the variance is infinite - such random walks 
are known as "Levy flights" - then the most efficient way to go far from the 
origin is to make one big step, accompanied by smaller steps. Note that the 
distance from the origin reached after a large number of steps has very different 
distributions in these two situations: Gaussian in the first opposed to a 

power-law tail in the second case. The interested reader may see [SJ, PP- 42-59, 
for a pedagogical introduction to Levy statistics. 

In the problem of independent multiple scatterings that we have discussed in 
this paper, the "step size" is the transverse momentum acquired by the gluon at 
each scattering, which has a probability distribution that falls like a{l±) ^ IJ^ 
in the MV model (even slower if there is an anomalous dimension different from 
unity). The variance of the step sizes, J cPl± l'± o-{l±), is thus infinite, and 
our problem falls in the category of Levy flights. Many of our results can be 
understood from this analogy. 
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